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Abstract. 

A kinetic theory of the Bohm criterion is developed that is based on positive- 
exponent velocity moments of the plasma kinetic equation. This result is contrasted 
with the conventional kinetic Bohm criterion that is based on a v moment of the 
Vlasov equation. The salient difference between the two results is that low velocity 
particles dominate in the conventional theory, but are essentially unimportant in the 
new theory. It is shown that the derivation of the conventional kinetic Bohm criterion 
is flawed. Low velocity particles can cause unphysical divergences in the conventional 
theory. These divergent contributions are avoided with this new approach. The two 
theories are compared using example distribution functions from previous presheath 
models. The importance of ion-ion and electron-electron collisions to determining the 
particle distribution functions throughout the presheath is also discussed. A kinetic 
equation that accounts for wave-particle scattering by convective instabilities is used 
to show that ion-acoustic instabilities in the presheath of low temperature plasmas 
(where T e 3> T<) can cause both ions and electrons to obtain Maxwellian distribution 
functions near the sheath. 
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1. Introduction 

The Bohm criterion, as originally derived pQ, provides a condition that the ion flow 
speed must satisfy as ions leave a quasineutral plasma and enter a nonneutral sheath. 
Bohm's derivation assumed a single species of monoenergetic ions with velocity Vj and 
Boltzmann electrons with temperature T e . From these assumptions, he showed that 
ions must be supersonic as they leave the plasma: Vi > c s = ^jT e /Mi pQ. It has since 
been shown theoretically [2] and experimentally [3] that equality typically holds in the 
Bohm criterion. Although the assumptions that Bohm made in his seminal work often 
provide an accurate approximation of physical systems, an important question to answer 
is: How does the Bohm criterion change for more general electron and ion distribution 
functions? 

Attempts to answer this question have been the topic of several theory papers over 



the past fifty years BEl IS1 IS1 IZl IH1 IS1 dOl HH 1121 HS1 031 H51 IISl HZl HSl HSl 1201 I2D 1221 



1231 [2H [22]- Harrison and Thompson put forth one of the first such generalizations of 
the Bohm criterion [4j. Like Bohm, they assumed that the electron density obeys a 
Boltzmann relation, but allowed for a general ion distribution, /j(v). They predicted 
that ions should obey the criterion J d 3 vv~ 2 fi(y)/ni < cj 2 , at the sheath edge. Later 
work by Riemann JTUJ [16] generalized the electron term as well and found the criterion 



Equation ([I]) is commonly called the "generalized Bohm criterion" or the "kinetic Bohm 
criterion." It is a frequently cited result that has been applied in both analytic and 
numerical models of laboratory [HI H21 031 Q31 HS1 001 HZ1 HH1 HOI 1201 1211 1221 1231 1231 
[25] 126] [27] and fusion [281 ES] plasmas. It is also discussed in a popular textbook [30J . 
In deriving equation (jl|), the electric field of the sheath is taken to be aligned in the z 
direction, and it is assumed that the only spatial gradients of the particle distribution 
functions are caused by this electric field. 

Although it is frequently cited in the theoretical literature, equation ([!]) does not 
produce a meaningful criterion for many plasmas of interest. If the ion distribution 
function has any particles with zero velocity, the left side of equation ([I]) diverges. If the 
velocity gradient of the electron distribution does not vanish at v z = 0, the right side of 
equation diverges. Consider the common example where the ion distribution function 
is a flow-shifted Maxwellian and the electron distribution is a stationary Maxwellian. 
In this case, the left side of equation is infinite and the right side is n e /T e , which 
produces an unphysical prediction. Equation places undue importance on the part 
of the distribution functions where particles are slow. 

These shortcomings of equation ([I]) have been pointed out before. Shortly after the 
publication of Harrison and Thompson's result, Hall discussed the fact that unphysical 
divergences arise when the ion distribution contains slow particles [5]. Lieberman and 
Lichtenberg (30] mentioned that equation (JT| can lead to "mathematical difficulties" at 
low energies. However, despite the fact that it often gives unphysical results and that 




(1) 
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this shortcoming has been discussed in previous publications, equation pi) continues to 
be used in plasma physics literature. 

Experimental and numerical studies of presheath physics have been published that 
appear to contradict equation 0. Laser induced fluorescence has been used to measure 
the ion distribution function in the presheath of laboratory plasmas [3TJ E21 E2] • These 
experiments found a Maxwellian ion distribution with a flow speed approaching the 
sound speed at the sheath edge; thus, fi(v z = 0) ^ 0. Numerical simulations by 
Sheridan [M], Meige et al [35] and Robertson [36] also find ion distribution functions 
with populations of low energy ions throughout the presheath. These used hybrid (fluid 
electron, kinetic ion) and fully kinetic particle-in-cell techniques. Consistent with the 
experiments, these simulations find that the ion flow speed approaches the sound speed 
at the sheath edge. These studies apparently contradict the implication of equation (HI) 
that the low energy part of the ion distribution function dominates the ion dynamics at 
the sheath edge. 

In section [2j we reconsider previous derivations of equation (JT1). We show that 
these derivations contain two errors. The first of these is taking the v~ l moment of the 
collisionless kinetic equation (i.e. Vlasov equation). Neglecting the collision operator is 
a mistake because the v ~ 1 moment of it diverges if the distribution functions contain 
particles near zero velocity. Not all distribution functions of interest satisfy the Vlasov 
equation. Only velocity moments with a positive power should be applied, otherwise 
divergences can arise in the integrand at v z — for non- Vlasov solutions. The second 
error is a mathematical mistake where integration by parts is misapplied to a function 
that is not continuously differentiable. This error can easily be corrected, but the 
resultant criterion then differs from equation (fij). 

In section |3j we derive an alternative kinetic Bohm criterion that is based upon 
positive exponent velocity moments of the kinetic equation, rather than the v ~ 1 moment 
of previous work. This approach avoids the possibility of divergent velocity-space 
integrals. Particles with low energy do not have any special significance in this theory. 
It also supports previous derivations of the Bohm criterion based on fluid equations, 
and it returns these results in the fluid limit. In contrast, equation ([I]) cannot return 
fluid results because it places undue importance on low energy particles. In section |4| 
we consider some example distribution functions that are common in low temperature 
plasmas, but for which the generalized Bohm criterion we derive in section [3] can give 
significantly different predictions than equation ([!]). We consider example presheath 
models with a cold ionization source (the Tonks-Langmuir problem), a warm ionization 
source as well as a collisional example where ions and electrons are Maxwellian. In 
section |5j we discuss how electron-electron and ion-ion collisions in the presheath are 
often important in order to determine the distribution functions near the boundaries of 
laboratory plasmas. Specifically, we show how ion-acoustic instabilities can significantly 
enhance these collisions and drive the distribution functions toward Maxwellians. 
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2. The conventional kinetic Bohm criterion 

2.1. The sheath edge and sheath criterion 

To derive equation (jll), one must first develop a mathematical definition of where the 
quasineutral plasma (presheath) ends and the nonneutral sheath begins. If the Debye 
length is shorter than the ion momentum transfer collision length (Ad Aj m ), the 
presheath and sheath are separated by a transition region [TB]. This transition region 
scales as A^A^ 5 [21 E], which is typically much shorter than the presheath length 
scale (I). In the asymptotic limit — > 0, the transition region is so short that the 
boundary layer problem can effectively be split into the two regions of quasineutral 
plasma (presheath and bulk plasma) and nonneutral sheath [161 El EH]- The sharp 
boundary that results in this limit is what is typically referred to as the "sheath edge." 

The sheath criterion is a condition that the spatial gradient of the charge density 
must satisfy at the sheath edge. Riemann derives the sheath criterion by applying to 
Poisson's equation the physical requirement that quasineutrality breaks down at the 
sheath edge [17] . His derivation can be summarized as follows. Expanding Poisson's 
equation about the sheath edge, where the electrostatic potential is referenced to zero 
= 0, yields V 2 = — 47r[p(0 = 0) + dp/d(p\^ =0 (j) + ...]. Here p = J2 S q s n s is the 
charge density. It is assumed that the electrostatic potential variation of the presheath 
and sheath is one- dimensional and we align this in the z direction. At the sheath 
edge, which is the boundary of quasineutral plasma, the first nonvanishing term in this 
expansion is the linear term: d 2 (f)/dz 2 = —Andp/d^l^o^- Multiplying this by d(f)/dz 
and integrating with respect to z yields 



E 2 dp 

h -jj 

47r d6 



i2 



C (2) 

0=0 



in which C is a constant. Since (p — > as z/Xn — > oo on the sheath length scale, the 
constant C must be zero [TTj . We are then left with dp/dipl^o = — E 2 / (^TKp 2 ) , which 
implies the sheath criterion 

1 £0 - < 3 > 

d(p 4>=o 

The sheath criterion provides a mathematical definition of the sheath edge. Using 
the fact that dn s /d(f) = —E~ l dn s /dz, the sheath criterion can also be written 
q s dn s /dz\ z= Q > 0. Since the relation between density and the distribution function 
is simply n s = J d 3 v f s , equation ^ becomes 

2.2. Derivation of the conventional kinetic Bohm criterion 

Equation ([T]) has been derived using various approaches by different authors [HUTT-I, but 
each approach contains the following elements. These are collisionless theories based on 
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the one-dimensional steady-state Vlasov equation 

V ^I± + 3^ E ^Il = o (5) 

z dz m s dv z 

where it is assumed that the only spatial gradients of f s are due to the electric field in 
the presheath and sheath. Each approach (at some point in the derivation) also divides 
equation ^ by v z , to obtain an expression for df s /dz, which can be put into the sheath 
criterion of equation Q. In other words, the v~ l moment of the Vlasov equation is 
used. Applying this approach yields the following form of a kinetic Bohm criterion 

E ^fV„I^<o. (6) 

^ m s v z dv z 

Assuming that the plasma consists of a single species of ions with charge qi = e and 
electrons, this is 



d 3 v _^l<__ / tfv-^l. (7) 

Mi v z dv z m e v z dv z 

One additional step is required in order to write equation ^ in the form of 
equation Q. This is to integrate the v z component of the ion term by parts, 



dv g —^= [ dv z ^-(— )+ [ dv z h- (8) 
v z ov z dv z \v z J v{ 

where the surface term (the first on the right side) is taken to vanish since [fi/v z ] Vz= ± oc = 
[U IB]. Applying equation QSJ) to the left side of equation (jT]) completes this simple 
derivation of the conventional kinetic Bohm criterion from equation Q. 



2. 3. Deficiencies of the conventional kinetic Bohm criterion 

Two mistakes are made in all derivations of the conventional kinetic Bohm criterion, 
including the one presented in section |2.2| They are: 

(i) The collision operator should not be neglected if one is to take the v~ x moment of 
the kinetic equation. This is because J d 3 v C(f s )/v z diverges unless the collision 
operator vanishes. However, the collision operator only vanishes when the plasma 
is in equilibrium, which implies that both ions and electrons have Maxwellian 
distributions with the same temperature and flow speed %F7\ EE1 EH1 HO]. This 
would lead to a contradiction in equation ([T]) because it reduces to oo < n e /T e for 
Maxwellian distributions. 

(ii) Since the function (l/v z )df s /dv z is not generally continuously differentiable, the 
integration by parts conducted in equation ^ is invalid. This deficiency has also 
been pointed out before [2T] . 

The easier of these two issues to correct is (ii) since the integration by parts step 
can simply be avoided and the kinetic Bohm criterion left in the form of equation ^ 
or ([7|. However, even these equations are incorrect because of issue (i), which will be 
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discussed below. That the integration by parts step of equation (|8J) is incorrect can be 
shown from a simple example. The contentious step is of the form 



dx— ( j-= [ dx-^-(—\ + / dx^, (9) 

X (XX J ^y-, (XX \X ' J X 



— 



for any physically possible distribution function (e.g., the restrictions /(±oo) = and 
that / is positive for all x can be imposed since any physical distribution must obey 
these). If one takes as an example / = exp(— x 2 ), the left side of equation fcty can be 
evaluated directly 



oo 1 jr poo 



oo 



dx- d t = -2l dxe~ x ' = -2^. (10) 
x dx 



However, if the surface term on the right side of equation (|9]) is taken to be zero, as it is 
in derivations of the conventional kinetic Bohm criterion, the right side of equation (|9]) 
diverges 



e ~ x2 ( Z" - ' 6 ' e~ x2 f°° e~ x2 



X 2 ^0\J^ X 2 ./[ e | x 2 



dx — — = limf / dx — r — h / dx — — ) (11) 

ar 

2 



- 2V7T + lim 



e" |e| +2v^erf(|e|) 
l lei 



7r + lim — e ' e ' 2 — > oo. 

The reason that integration by parts cannot be applied in equation (J8| is that it is 
only valid for continuously differentiable functions [43J. However, v~ l df/dv z is not 
continuous unless df/dv z \ Vz= o = 0, and v~ 1 df/dv z is not continuously differentiable 
unless both v~ 1 d 2 f/dv 2 and v~ 2 df /dv z are continuous. Many physically reasonable, 
often expected, plasma distribution functions do not satisfy these properties. Thus, issue 
(ii) restricts the previous kinetic Bohm criteria to the form of equation (JTl). However, 
issue (i) will show that there is a problem with equation (J7|) as well. 

Equation (Q also places undue importance on the low-velocity part of the 
distribution functions. The primary deficiency of the collisionless Vlasov approach is 
simply that the collision operator cannot be neglected if one is interested in v~ x moments 
of the kinetic equation. To illustrate this, consider what happens if the approach of 



section 2.2 is taken but the collision operator is not neglected. Then, the relevant 



kinetic equation has the form 

v*7r + —E?fi- = C<J.), (12) 
oz m s ov z 

in which C(f s ) = S S 'C(/s;/s') is the total collision operator. The total collision 
operator consists of the sum of component collision operators describing collisions 
between the test species s with all species in the plasma (s f ) including itself (s = s'). 
The component collision operators have the Landau form 



C(f s ,f s/ ) = ~ ■ f dV Q- ~ —4-)f.{v)f*W, (13) 

(7v J \m s i ov m s o\y 
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in which Q is a tensor kernel. Lenard [3?] and Balescu [42] have shown that in a stable 
plasma Q = Qlb where 

kk5[k- (v- v')] 



Q 



LB 



2« 



tffc- 



fc 4 |efk,k-v) 



Here 



^ k 2 m s J 



d 3 v 



k ■ df s /dw 
kv 



(14) 



(15) 



is the plasma dielectric function for electrostatic fluctuations in an unmagnetized plasma. 

Taking the v~ x moment of equation (12), in order to find an equation for df s /dz, 
and putting the result into equation Q gives the criterion 



E 



a v- 

-oo 



v z dv z 



< 



E 



E 



OO 1 

'3 1 



d'v-Cifs 

V z 



(16) 



which can be compared to the Vlasov result from equation (|6]). A brief study of 

ide, whicl 

.1^=0 ^ or f s (v z 



equation (16) shows that not only the left side, but also the right side, which depends on 

0)^0 



the v~ l moment of the collision operator, diverges if df s /dv z 
for any species s. 

Equation (16) shows that neglecting the collision operator is not a consistent 
approximation when the v ~ 1 moment is applied. For example, consider a plasma with 
a stationary Maxwellian electron species and a Maxwellian ion species flowing relative 
to the electrons. In this case C(fi, fi) = and C(f e , f e ) = 0, but C(/ e ,/i) ^ 0. Since 
fi(v z = 0) 7^ 0, the C(/e, fi) term will cause the right side of equation (16) to diverge. 
The ion term on the left side of equation (16) diverges for this example as well. The 
only case for which the Vlasov approach, and thus equation ([7]), is strictly correct is 
for Maxwellian ions and electrons with the same temperature and flow speeds: In this 
limit C(f e ,fi) = 0. For this case, equation ^ [or equivalently equation (16)] reduces 
to n/T — n/T < 0, which is a true statement, but is not useful as a Bohm criterion. 

The presence of an ionization source in the Vlasov equation can produce additional 
problems for the conventional kinetic Bohm criterion. Including the effects of a source 
Si in the above analysis, one finds an equivalent divergent contribution to the right 
side of equation (16) unless Si(v z = 0) = 0. Section 4.2.1 discusses the implication for 
the Bohm criterion of various ionization source functions that have been considered in 
previous literature. In some of these examples, such as a Maxwellian source function, 
the resultant ion distribution function contains slow particles causing the conventional 
kinetic Bohm criterion to diverge. 



3. A kinetic Bohm criterion from velocity moments of the kinetic equation 

By taking positive velocity moments (e.g., v, v 2 , vt> 2 , . . .) of the full kinetic equation, 
a formally exact set of fluid equations are derived. However, this approach suffers from 
the "closure" problem where evolution equations for low order moments are determined 
by higher order moments. Because of this accounting difficulty, the fluid equations are 
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not closed unless further approximations are made. In the limit that the distribution 
function is a flow shifted Maxwellian, the plasma is completely described by the lowest 
order moments (1, v, and v 2 ), which determine the evolution of n s , V s and T s . However, 
deviations from Maxwellian can require evaluation of higher order moments in order to 
track the evolution of each parameter specifying the distribution. 



3.1. Fluid moments of the kinetic equation 

A hierarchy of fluid moment equations can be constructed from velocity moments of the 
plasma kinetic equation for species s 

f + v.^ + ^E.^ = C(/.) (17) 
ot ax m s c>v 

by applying the standard definitions of the fluid variables in terms of velocity-space 
integrals of the distribution function. These are density 



oo 



oo 



n s = I d v f a 
fluid flow velocity 

V, = — / d 3 v vf s , (19) 



n 

scalar pressure 

'■OO 



stress tensor 



oc 



oo 



temperature 



p s = I d 3 v -m s vl f s = n s T s , (20) 

J-oo & 

U s = J d 3 vm s (v r v r -^v 2 r l^)f S} (21) 



1 f°° . 1 1 

T s = — I d 3 v -m s v r f s = -m s v Ts , (22) 
n. . n 3 2 



and frictional force density 

poo 

R s = / d 3 v m s vC{f s ). (23) 



Here we have defined the relative velocity v r = v — V s , where V s is the fluid flow 
velocity from equation (fl9j). 



The density moment (J d 3 v . . .) of the kinetic equation (17) yields the continuity 
equation 

£ + ».(«.V.)-0. (24, 
The momentum moment (J d 3 vm s \ . . .) yields the momentum evolution equation 

m s n s + V s • — - = n s q s E - - — • n, + R s . 25 

V dt dx J H 9x 9x v ' 
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Equations (24), (25) and subsequent equations built from higher exponent velocity 
moments of the kinetic equation constitute a hierarchy of fluid equations. In the next 
section, we use equations (24) and (25) to formulate a Bohm criterion that can be 
written in terms of the distribution functions by associating the fluid variables with 
their definitions in terms of f s from equations (18) - (23). In this way the resultant 
criterion retains a kinetic form. 



3.2. A kinetic Bohm criterion 



Bohm's original criterion (V* > c s ) was a condition concerning the ion speed (assumed 
to be monoenergetic in his paper) at the sheath edge. Kinetic Bohm criteria seek to 
generalize this condition to account for arbitrary ion and electron distribution functions. 
However, it is unclear exactly what qualifies as a "kinetic Bohm criterion." For instance, 
the sheath criterion of equation Q specifies a condition that the spatial gradients of 
the distribution functions must satisfy at the sheath edge, yet it is not typically called 
a "Bohm criterion." A possible definition for "Bohm criteria" might be statements 
concerning the ion flow speed at the sheath edge. Although it is not obvious that the 
conventional kinetic Bohm criterion of equation always fits this definition, we will 
look specifically for a condition concerning the fluid flow velocity of ions [defined in 



terms of equation (19)] at the sheath edge. 



We assume that the plasma is in steady-state and that the only spatial variation 
in f s is due to the electric field in the sheath and presheath. We take this electric field 



to be in the z direction. 



Thus, the fluid variables n s , V s , T s , p s , U s and R s are only 



functions of the spatial variable z. With these assumptions, the continuity equation (24) 



and momentum equation (25) reduce to 

dV z . 



dz 



+ V Z . 



dn s 
' dz 







and 



m s n s V ZyS 



dV z . 



n s q s E - 



dp s 



dU, 



+ R 



(26) 



(27) 



dz ' ' ' dz dz 
in which the z subscript refers to the z component of a vector and the subscript zz to 
the zz component of a tensor. 



Solving equation (26) for dn s /dz and putting the result into the sheath criterion 
J2 S q s dn s /dz\ z=0 > [equation (HI) 



dV z . 



' V z , dz 



yields 
< 0. 



2 = 



(28) 



Equation (28) is a condition concerning the spatial gradient of V s at the sheath edge. We 



are interested in a condition on V s itself, which can be obtained by using equation (27) 



to find an expression for dV ZjS /dz. Putting equation (27) into (28) yields 

q s n s - (n s dT s /dz + dH ZZjS /dz - R Z)S )/E- 



m s Vl s - T s 



< 0. 



J 2 = 



(29) 
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Equation (29) is a kinetic Bohm criterion that provides a condition that the flow 
speed of ions must satisfy at the sheath edge. It makes no assumptions about the 
distribution functions and can be written explicitly in terms of them by substituting 
the fluid variables with their definitions from equations (18) - (23). It does depend 
on spatial gradients of the higher-order moments temperature and stress tensor and 
the collisional friction. These could be eliminated in terms of spatial gradients of even 
higher order fluid moments, the heat flux in this case, by using the temperature evolution 
equation [obtained from the (J d 3 vv 2 . . .) moment of the kinetic equation]. However, no 
matter how far one carries out the hierarchy expansion, the subsequent Bohm criterion 
will still depend on a spatial derivative of f s inside of some fluid moment integral. 
Various closure schemes can be applied to deal with these higher-order moments [23] . 

In many plasmas, but not all (see e.g. section 4.2.2), the temperature and stress 
moments vary on spatial scales much longer than the Debye length. For instance, the 
gradient length scales for the temperature and stress terms are often characteristic of 
the length scale for collisions between species s and the background neutrals or another 
species (s' ^ s). These lengths are typically on the order of the presheath length scale 
(I). Equation (29) can be simplified significantly for such plasmas because the terms in 
parentheses become negligible. The friction term is also controlled by collisions and acts 
over a collision length scale (including friction with neutrals [44]). One the other hand, 
the gradient length scale of the electrostatic potential approaches the Debye length at 
the sheath edge. Thus, the terms in parentheses in equation (29) are 0(Xd/1) <C 1 
smaller than the q s n s term in such plasmas. In fact, many presheath models are based 
on the limit \d/1 —> 0, in which case the electric field is found to become infinite at the 
sheath edge and the terms in parentheses in equation (29) formally vanish [T21 EES]. If 
these terms are negligible, equation (29) reduces to 



E 



m sV Z 2 g - Tg 



< 0. 



(30) 



Equation (j30J) can be further simplified if the electron fluid flow speed is much 
slower than the electron thermal speed at the sheath edge (V Zje <C VT e )- In this common 
situation, equation pOl reduces to 



E 



< 1, 



(31) 



in which i label the different ion species. Equation (31) was first derived by Riemann 
using a fluid approach [17]. The limit of equation (31) is obtained from the general 
Bohm criterion [equation (29)] by effectively asserting that deviations from Maxwellian 
plasmas are small and the conventional fluid theory is valid. Equation (31) can be 
written explicitly in terms of f e and fa by substituting equations (18), (19) and (22) for 
n s , V s and T s . 

To connect with Bohm's seminal work, consider a plasma with a single species of 
ions of unit charge. Then equation (31 ) yields 

V z>i > Jcl + vy2. (32) 
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If T e ^> Tj, equation (32) simply reduces to the usual Bohm criterion V z ^ > c s . However, 
whereas Bohm assumed monoenergetic ions, here V Z) i and the T e in c s = ^jT e /Mi are 
defined in terms of velocity-space moments of the ion and electron distribution functions 



[equations (19) and (22)]. Explicitly in terms of f e and /j, the V Zj i > c s criterion can be 
written 



d 3 vv z fi > 



d 3 vf 



d 3 vv 2 J ( 



1/2 



(33) 



in which / in the j d 3 vf term can be either fa or f e because of quasineutrality at the 
sheath edge (n e ~ n^). 

4. Example distributions for comparing the different Bohm criteria 

In this section, we use ion and electron distribution functions from example presheath 
models discussed in previous literature to demonstrate the similarities and differences 
between the conventional kinetic Bohm criterion and the criterion derived in section [3] 
The Debye length is much shorter than the ion or electron momentum transfer collision 
length in all of these examples. They are also concerned with single ion and electron 
species plasmas. Thus, we will be comparing the conventional Bohm criterion of 
equation Q with the single ion and electron species version of equation (29) in 
which the densities, flow speeds and temperatures are calculated from and f e using 
equations (18), (19) and (22). The example ion distribution functions considered 



include a monoenergetic distribution, collisionless models including the Tonks-Langmuir 
problem and variations of it that account for finite temperature ionization sources, as 
well as a Coulomb collisional presheath model where ions have a flowing Maxwellian 
distribution. The example electron distributions considered are Maxwellian and a 
collisionless model where the Maxwellian is truncated for velocities corresponding to 
electrons that have escaped the plasma through the sheath potential drop. 



4-1. Monoenergetic ions and Maxwellian electrons 

The idealized plasma that Bohm considered in his original paper [TJ assumed 
monoenergetic ions, /, = nj<5(v — Vi), and Maxwellian electrons 



fe 



7?.,, 



3/2 ?; 3 



71 



Te 



U 



Te 



(34) 



For these distribution functions, the components of the conventional kinetic Bohm crite- 
rion of equation Q are / d 3 v fi/(MiV 2 z ) = ^/(M;^ 2 ) and -/ d 3 v (df e /dv z )/(MiV z ) = 
n e /T e . Putting these into equation Q and applying the quasineutrality assumption 
gives Bohm's original criterion Vi > c s . Equations (21) and (22) give Tj = and 
n^ 2j i = for monoenergetic ions. For Maxwellian electrons, V, = 0, Il 22je = and 
dT e /dz = (see section 4.2.2). Thus, equation (31) reduces to the same condition 



Vi > c s . The conventional kinetic Bohm criterion and the kinetic criterion developed 
in section [3] both reduce to the criterion of Bohm's original work [1] for monoenergetic 
ions and Maxwellian electrons. 
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4-2. Collisionless presheath models 

4-2.1. Collisionless ion effects Ion kinetic effects in the presheath are often studied 
using hybrid models in which ions are treated with a collisionless kinetic equation, 
and electrons are assumed to obey a Boltzmann density profile, n e = n o exp(e0/T e ) 
[H El HH H21 H31 [221 ESI EH ESI S5]. The main distinguishing feature in the various 
models is how the ion source term is treated. The first of these theories was proposed in 
the seminal work of Tonks and Langmuir [15] , later extended by Harrison and Thompson 
[I] . In this model ions are assumed to be born with zero energy; all of their energy being 
provided by subsequent acceleration through the presheath electric field. Riemann has 
shown that the ion velocity distribution function at the sheath edge in the Tonks- 
Langmuir problem [where the ion source is singular in velocity Si oc is independent 

of the ionization rate (22] • It satisfies the conventional Bohm criterion of equation ([!]) 
[U |22], which is possible because fi(v z = 0) = at the sheath edge. 

The Tonks-Langmuir problem has also been generalized to account for ion source 
distributions that have an energy spread. Emmert et al [9] assumed a source function 
of the flux form 

S { = h(z)^-exp(-v 2 Jv Ti ). (35) 

V Ti 



in which h(z) is the spatial strength of the source. Equation (35) vanishes for v z = 0, 
so the v ~ 1 moment does not diverge. A consequence of this source distribution is that 
the resultant ion velocity distribution function satisfies fi(v z < 0) = at the sheath 
edge [9]. Elsewhere in the presheath fi(v z < 0) ^ 0. Bissell has shown that Emmert's 
solution satisfies the conventional Bohm criterion of equation [12J, which is only 
possible because Si(v z = 0) = 0. 

Since ions are born from ionization of neutral atoms, which themselves are expected 
to have a Maxwellian distribution, a natural assumption might be that the ion source 
function also have a Maxwellian distribution with a temperature near the neutral gas 
temperature. Bissell and Johnson [HI [T3] applied an analysis similar to Emmert's, but 
assumed a Maxwellian 

^ = ^exp(-^/4,)- (36) 

VTi 

Using this Maxwellian source function, the resultant ion velocity distribution function 
at the sheath edge was found to have a nonzero contribution at v z = 0. Sheridan [31] 
has also provided a detailed numerical study of Bissell and Johnson's model showing 
that fiiy z = 0) 7^ at the sheath edge for any finite ion temperature in the source 
distribution. Robertson has found similar results for warmer source temperatures than 
Sheridan considered Similar ion velocity distribution functions with fi(v z = 0) ^ 
at the sheath edge have also been found in particle-in-cell simulations [24J, and other 
kinetic simulations [TJ]. Since fi(v z = 0) ^ at the sheath edge for this model, the 
left-hand side of equation ([TJ diverges and the conventional kinetic Bohm criterion does 
not provide a useful condition limiting the ion speed at the sheath edge. However, if 
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the Bohm criterion from section [3] is applied to these ion distributions, we find that the 



essentially fluid result from equation (32) holds, where the fluid variables are defined in 



terms of fi from equations (19) and (22). 



4-2.2. Collisionless electron effects If the electron momentum transfer collision length 
is longer than the sheath length, the electron distribution function at the sheath edge 
will be depleted for velocities exceeding some threshold. This threshold corresponds 
to the velocity required to overcome the sheath potential drop. Electrons that escape 
through the sheath are lost to the boundary surface, so the electron distribution function 
is depleted beyond this critical velocity by a factor of order d/X e , where d is the distance 
to the boundary and A e is the electron collision length. If the material boundary is in 
the +z facing direction, the truncation velocity is given by vii c = —^2e(\(j)b\ + 4>)/m e z. 
Here <f>\, is the potential of the boundary surface with respect to the plasma. We have 
chosen the plasma potential as the reference potential (0 P = 0), so <fi typically takes 
negative values through the presheath and sheath. 

Assuming that electrons are Maxwellian, aside from the truncation, the electron 
velocity distribution function can be written 

— 2 

f* = 3ifer ex pf-H2-)^(^ + w n, c ) (37) 



7T 



J Te x u Te ' 

in which H is the Heaviside step function, VTe 

of the fluid variable definitions from equations (|18|)-(j22|), the density is 



^/2T e /m e , and v\\ )C = |v[| )C |. In terms 



n. 



1 + erf f ^ 



the flow velocity is 
V e = 

the temperature is 

T„ = T <! 1 - 



1 exp(-ujj |C /t# e ) 
7r 1 + eTi(v\\ }C /v Te 



;V T e Z, 



eX P(- V lc/ V Te) 



2 exp(-2v| c /4J 
3a/tt v Te [1 + erf(t)|| )C /w Te )] 3vr [1 + erf(w|| jC /u Te )] ; 
and the zz component of the stress tensor is 



-n P T R 



expl — 



J Te 



V\\,c 
VTe 



+ 



1 exp(-i| e /4 e ) 
^/tx 1 + erf (dm c /w Te )J 



(3f 



(39) 



(40) 



(41) 



Equations (38 )-(41 ) are shown in figurenlas functions of v\\ tC /vT e . In this problem, we will 



use the temperature evolution equation to close the moment equations. Thus, the heat 



flux, q s = m s J d 3 v\ r v^f s /m s , will be required. For the distribution of equation (37) 
the electron heat flux is 

.2 T T T T2 - 



n e V e TJ -- + 



v 



V. 



II c V 

1^ + 3, 



V P 



+ 2^- 



Te 



VTe VTe 



J Te 



(42) 



Figure (fll) shows that if v\\ c /vT e ^ 1, the dT e /dz and dH ZZte /dz terms of 



equation (29) may become important. This is a consequence of the distribution being 
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0.5 1 1.5 2 2.5 



v \\,c/ v Te 



Figure 1. Normalized electron density (solid black line), flow speed (dashed red line), 
stress (dash-dotted green line) and temperature (dotted blue line) computed from the 
moment equations (18), (19), (21) and (22) for the truncated Maxwellian distribution 
function of equation (37). 



far from Maxwellian. In this circumstance, the simplification of equation (30) is not 
expected to be valid. To demonstrate that the temperature and stress gradients 
should not be neglected when the distribution is far from Maxwellian, we compare 
the predictions of equations (29) and (30). Assuming that ions are monoenergetic, 
equation (30) reduces to 

V? > c s ^l-m e V e 2 /T e . 



(43) 



Putting equations (38), (39) and (40) into (43) yields the following form of a Bohm 
criterion 



exp(-vjf A4J 8 exp( 



3^v Te [1 + erf(u|| iC /ur e )] 3vr [1 + erf(v\ Uc /v Te )f 



1/2 



.(44) 



Equation (44) is shown as the red dashed line in figure |2| It is an essentially fluid 
result which neglects the dT e /dz and dYi ZZ)e /dz terms that must be determined from 
higher-order moments of the kinetic equation. 

A full kinetic calculation of the Bohm criterion from equation (29) requires taking 
the spatial derivative of equations (40) and (41). At this point, both n e and T e in 
equation (37) are expected to have a spatial dependence. A closure scheme is required 
to solve for the three variables n e , V e and T e . This can be provided by deriving a 
temperature evolution equation from the J d?vw r v^ moment of the kinetic equation (12). 
The 1-D steady-state version for electrons, neglecting collision terms, is 
d_ 
dz 



Qz,e + 



^n e T e + -m e n e V* ( 



V z . e + V Z M 



+ en e V z , e E = 0. 



(45) 



From equations (38)-(42), we find that 



Qz,e + 



^n e T e + -m e n e V* e 



Vz.e + V Z M 



z,e ±J -zz,e 



n e V e T e [2 + 



VTe 



(46) 
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v \\,c/vTe 



Figure 2. Prediction for the minimum ion speed at the sheath edge normalized 
to c s as a function of the truncation speed for the electron distribution function of 
equation (37). The red dashed line shows the prediction of equation (44), which is an 
approximation that assumes electron temperature and stress gradients are negligible. 
The black solid line shows the prediction of equation (53), which is the full kinetic 
result. 



Putting this into the collisionless temperature evolution equation (45), and identifying 
the continuity equation d(n e V Zje ) / dz = 0, yields 

„,2 



TJ2 + 



Noting that 



d_ 

dz . 



dz \ Vj, e 



-eE. 



J Te' J 



d 

dz 



e(\<p b \ + <f>) 



eE 



1 dT e 



v^ e T e dz 



equation (47) reduces to 
df e 
dz 



0. 



(47) 



(4£ 



(49) 



Thus, in the absence of collisions, the parameter T e is constant. The temperature T e , of 
course, changes in space as vu c varies. However, it is also useful for the other examples 
to note that for a full Maxwellian T e — > T e , so equation (49) shows that the temperature 
is constant for these examples when electron collisions can be neglected. 



Equation (49) implies 

dz \VTe' 



1 VTe eE 



(50) 



Using equation (50) to evaluate n e dT e /dz + dH ZZfi /dz from equations (40) and (41), one 
can show 



1 / dT P dll, 



en e 



E 

en ( 



n e -r- + 
dz 



dz 



(51) 



(m e V e 2 -T e ) 1 + 



VTe 



exp(-vlJv^ e ) 



U|| )C V^[1 + erf(« ||, c /t>Te)] 
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Putting this into equation (29), along with the monoenergetic ion term, yields 

exp(- 



e 2 nj 

M;V? 



1 + 



VTe 



In 



7r[H-erf(uii c /«r e )] 



< 0. 



Rearranging, the full kinetic Bohm criterion of equation (29) yields 

,2 



K kin > C f 



1 + 



VTe 



e M-vlJv Te ) ,-1/2 



a/7t[1 +erf(ui| iC /ur e )] 



(52) 



(53) 



< 1 



equations (44) and (53) give 



Equation (53) is shown as the black line in figure [2 

As expected, figure [2] shows that for v\\ tC /vTe 
significantly different predictions. This is because the electron distribution function 
deviates significantly from Maxwellian in this region, and the dT e /dz and dH ZZte /dz 
terms become important. For v\\ c /vTe ^ 1, the conventional result that ions obtain 
the sound speed at the sheath edge is returned. For a floating boundary, the ion 
and electron fluxes balance at the sheath edge, and the sheath potential drop is 
\A<fi b \ = T e [l + ln(Mj/27rm e )]/2e. This is typically a few electron temperatures, in 
which case v» c /vTe > 1- Thus, the effects of a truncated electron distribution do not 
significantly affect the Bohm criterion near a floating boundary. 

A simple way to construct a boundary for which v» c /vT e ^ 1 is to bias a probe 
near, or more positive than, the plasma potential. Depending on the bias, and size of 
the probe relative to the chamber wall, the sheath near the probe can be an ion sheath, 
double sheath or electron sheath |46j. The parameter v\\ tC /vT e can be varied using the 
applied potential. Figure [2] shows that temperature and stress gradients are important 
in this regime and that one must use the full kinetic result from equations (29), rather 
than the common approximation from equation (30). 

Next, we compare the Bohm criterion of equation ( [53] ) to that predicted by 
the conventional kinetic Bohm criterion of equation Q. The electron term in the 
conventional criterion is 



1 



Idfe 



1 + 



exp(-t| c A4 e ) VTe 
?[l + erf(u|| iC /v Te )]u||, c 



Applying equation (54) and that ions are monoenergetic to equation ([I]) yields 

V Te eX P(-*|c/4e) 1 ~ 1 / 2 



Vi>c s il + 



7r[l + eri(v\\ )C /v Te )] 



(54) 



(55) 



This is the same as equation (53), which was obtained from the method of positive- 
velocity moments. 

In this example, the conventional kinetic Bohm criterion and the criterion from 
section [3] gave the same result because the electron distribution chosen is a possible 



solution of the Vlasov equation (it satisfies df e /dv z 



0). However, if the 



distribution function of equation (37) is modified slightly, the conventional result can 
change dramatically. For example, if the electrons have a small drift velocity in the z 
direction [so that v 2 — > (v — Y z >e ) 2 in the exponential of equation (37)], then equation M 
leads to the prediction l/(M i V i 2 ) < — oo. With the new method from section adding 



Kinetic Bohm Criterion 



17 



a small drift leads to 0(V e /vre) corrections to equation (53). These are negligible as 



long as the applied drift is much slower than the electron thermal speed. 
4-3. Collisional presheath 

Finally, we consider plasma in which the ion-ion and electron-electron collision lengths 
(A*~* and \ e ~ e ) are much shorter than the presheath length (which is typically the ion- 
neutral collision length). In this situation, the ion distribution function at the sheath 
edge is a flow-shifted Maxwellian. If the electron-electron collision length is longer 
than the sheath thickness, the electron distribution function at the sheath edge will 
be truncated at a velocity corresponding to the sheath energy. However, we found in 



section 4.2.2 that this truncation rarely affects the Bohm criterion. Thus, we assume 
that electrons have a stationary Maxwellian distribution. 

For flowing Maxwellian ions and stationary Maxwellian electrons, the Bohm 



criterion from equation (30) simply reduces to V; > \J c 2 s + v\ i j'l. However, the ion 
term in the conventional kinetic Bohm criterion from equation diverges for a flowing 
Maxwellian distribution 



[see equation (11)]. Thus, equation ([TJ leads to an incorrect statement, 00 < n e /T ei 
whereas equation p0| ) gives a criterion identical to that provided by conventional fluid 
theory. In section we show that collisional presheaths (A l_i ,A e_e <C I) are common in 



laboratory plasmas if the electrons are much hotter than ions. 
5. Coulomb collisions in the presheath 

In this section, the role of Coulomb collisions in the presheath is discussed. We are 
interested in determining when the ion-ion collision length is shorter than the presheath 
length. If this condition is met, the ion distribution function is expected to be a flowing 
Maxwellian at the sheath edge. This situation provides an example that is commonly 
found in laboratory plasmas, but where the conventional Bohm criterion and the one 



derived in section 3.2 provide very different predictions. 



Coulomb collisions in a plasma are described by the collision operator of 



equation (13). If the plasma is stable, Lenard [37] and Balescu [42] showed that the 



appropriate collisional kernel, Q, is equation (14). Recently, we generalized Lenard- 
Balescu theory to also account for unstable plasmas [391 HO] • For an unstable plasma, 
the collisional kernel is the sum of the stable plasma (Lenard-Balescu) term and an 
instability-enhanced collision term: Q= Qlb + Qie- The instability-enhanced collisional 
kernel is [391 SD] 



m. s ./ /c 4 ^ (lorj - k • v) 2 + 7| 
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exp(2 7i t) 

[(u Rj -k.vy + ^}\di(k } co)/dco\i; 

in which urj is the real part of the angular frequency and jj the imaginary part of the 
angular frequency of the j th unstable mode. Time in equation (57) is calculated in the 
rest frame of the unstable mode. For convective instabilities, t can be translated into 
distance in the laboratory frame [39]. We will be considering ion-acoustic instabilities 
here, which are convective. 

It has been shown that both the Lenard-Balescu J3JJ EH] and instability- 
enhanced [ID] terms cause component collision operators to evolve to uniquely 
Maxwellian distribution functions. Thus, the distance over which a non-Maxwellian 
distribution function evolves to Maxwellian is determined by the shorter of the stable 
plasma collision length A^q 8 ~ v/v^ s or the instability-enhanced collision length 
Ajg S « v/v^ s . Here v is the collision frequency, s is the particular species considered, 
s — s refers to self collisions within this species, and v is the velocity of the test particle 
considered in the rest frame of the s distribution (i.e., it is how far in velocity space 
a test particle must go to be thermalized, which is approximately |v — V s |). We are 
interested in like-particle collisions i — i and e — e, but not unlike-particle collisions e — i 
or i — e because they are much less frequent. 

The collision frequency can be estimated directly from the collision operator: 
u s ~ s ~ C(f s ,f s )/f s . Using d/dv ~ v" 1 in equation (13) shows that the collision 
frequency also consists of the sum of stable plasma and instability-enhanced terms 

C(f s , f 8 ) n s 

rvj 

fa m sV 2 

We first consider the stable plasma contribution, ~ n sQui / ( m s v 2 ), which 
requires evaluating Q^ s from equation (14). Since the dielectric function at u = k- v is 
adiabatic over most of the fc-space integral, e(k, k- v) ~ 1 + k^X^, the Lenard-Balescu 
kernel reduces to the Landau collisional kernel [H] 

m s u 6 

in which u = v — v' and A w 127rn e A|, e . Taking u ~ v, the dominant term of the stable 
plasma collisional kernel for like-particle collisions is Q^ 1 ~ 27rgjln A/(jn s v), and the 
stable plasma collision frequency is 

^LB S -^lnA. (60) 

Next, we consider the possibility of ion-acoustic instabilities enhancing collisions 
in the presheath. Ion-acoustic waves have a phase speed that is slow compared to 
the electron thermal speed, u/kvT e 1 and fast compared to the ion thermal speed 
(u — k- Vi)/kvTi ^> 1. We assume that electrons are much hotter than ions, Tj/T e 1, 
in which case ion Landau damping can be neglected. Under these assumptions, the 
electrostatic dielectric function of equation ( 15 ) reduces to 

1 U pi . V 71 " U 

£=1+ PAL " (W - k • V,)2 + VAL kVTe' (61) 



s-s w \Jsi JsJ s ( ns-s , n s ~ s \ /ct 

V 1 ~^l^LB + ^IE )• l 5 ^ 
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Figure 3. Real and imaginary parts of the ion-acoustic dispersion relation, 
equation (62), for three values of the ion fluid speed (normalized to the sound speed) 
that are found in presheaths: 1/4 (green, dash-dotted line), 1/2 (blue, dashed line) 
and 1 (red, solid line). We have used k • Vj ~ kVi in these plots. 



In obtaining equation (61), we have assumed that the electron distribution is 
Maxwellian. Accounting for a truncated Maxwellian of the form of equation (37) 
leads to negligibly small corrections as long as v\\ tC /vT e ^ 1- A dielectric 
function for the truncated Maxwellian can be developed using the incomplete plasma 



dispersion function jl?], which shows that corrections to equation (61) are of order 
0{exp(— v m c /vj, e )vT e /v\\,c}- We assume these are negligible here. 



Solving for the roots of equation (|61|) gives the ion-acoustic dispersion relation 



0J± 



k • V* ± 




\ 2 



Uj ^Jirm e /8Mi \ 
~ l + k 2 X 2 De fl 2 )- 




(62) 



y/l + k 2 X 

A growing wave is present as long as the ion fluid speed is large enough: |k • Vj 

3] for three representative 



> 



kc 



[1 + k 2 X 2 De ). Equation (62) is plotted in figure 



asma is neutral (m 



values of the ion fluid speed in the presheath assuming the p 
Figure [3] shows that the relevant wavelengths of unstable modes are near the electron 
Debye length (or shorter). 

Appendix A provides a calculation of the instability-enhanced collisional kernel 



when ion-acoustic instabilities are present. Putting this result from equation (A. 13) 



into equation (58), we find that the instability-enhanced collision frequency for like- 



particle collisions is 



'LB 



1+K 



IE 



81n A (1 + k 2 ) 2 



exp 



irm e ni Z 
16Min e X De 



when ion-acoustic instabilities are present. In equation (63) 



LB 



(63) 

is the stable plasma 



collision frequency from equation (60) and k c , defined in equation (A. 10), is a parameter 



that accounts for the fraction of fc-space that is unstable. Here Z is a spatial coordinate 
alined with z, but is shifted so that its origin corresponds to the location where the first 
unstable wave is excited. For the presheath, Z = is at the presheath-plasma interface, 
and Z takes positive values through the presheath. 



Finally, we apply equations (60) and (63) to the parameters of an experiment 
in which Claire et al [33] used laser induced fluorescence to diagnose the ion velocity 
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distribution function throughout the presheath. Similar experiments have also been 
conducted by others [3H [32]. Claire et al studied an argon discharge with plasma 
parameters T, = 0.027 eV, T e = 2.5 eV and n = 5.5 x 10 9 cm -3 at a neutral 
pressure of 1.8 x 10~ 4 mbar. They found that the ion velocity distribution function 
was Maxwellian in the entrance (bulk plasma facing) region of the presheath, had a 
three-part distribution in the middle of the presheath that was qualitatively similar to 
collisionless presheath models such as those discussed in section 4.2 (they compared to 
Emmert's [9] model in particular), and became Maxwellian again near the sheath edge. 
The finding that the ion velocity distribution function became Maxwellian between the 
mid presheath and the sheath edge was purported to be a new and surprising result that 
is not predicted by previous kinetic theories of the presheath. Here we show that this 
measurement can be explained by ion-ion collisions enhanced by convective ion-acoustic 
instabilities. 

Consider the stable plasma contribution to the ion-ion collision length, A^ ~ 
u/z/lb\ using equation (60). For the parameters of this experiment, A^ = 1.8 x 
10~ 13 v A m, in which v is in m/s. Near the entrance to the presheath, the ion flow 
speed is very small and the relevant speed is the ion thermal speed v ~ vt%- Putting 
in v = vxi = 3.4 x 10 2 m/s gives X££ = 2.4 mm. Near the presheath-sheath boundary, 
the ion flow speed is nearly the sound speed, which is much larger than the thermal 
speed, so v ~ c s . Using v = c s = 2.4 x 10 3 m/s gives A^ = 6.0 m. The presheath in 
these argon discharges is a couple of times the ion-neutral collision length [3] 08] . We 
consider the presheath to be collisional for ion-ion collisions if the ion-ion collision length 
is smaller than the ion-neutral collision length: /X~ n < 1. If this is the case, ion-ion 
collisions are frequent in the presheath and the ion distribution function should be a 
flow-shifted Maxwellian. At 1.8 x 10 -4 mbar, the neutral density is n n = 4.3 x 10 12 cm -3 . 
The total ion-neutral collision cross section for the energies of interest is approximately 
a ~ 1 x 10~ 14 cm 2 [30]. Thus, the ion-neutral collision length is A l ~ n ~ 1/ (n n o~) = 23 cm. 
According to the stable plasma contribution to the collision frequency, ions near the 
entrance to the presheath are collisional, but ions near the presheath-sheath boundary 
are not. 

An estimate of ion collisionality throughout the presheath can be made using a 
model for the ion flow speed. For this estimate, we apply a well known modified- 
mobility-limited-flow model from Riemann [2] 

e A=\J^L) and dz = dV i ( ( l a ~ V * \. (64) 

Two models for the ion-neutral collision length are commonly used: constant collision 
length \ l ~ n = I, v l ~ n = Vi/l, or constant collision frequency v % ~ n m c s /l. For simplicity, 
we apply the constant collision frequency model. The model of equation (64) has been 
verified experimentally to a distance z ~ 2A l_n jHlHS]- We choose a coordinate system 
where z = is the sheath edge, and z takes negative values throughout the presheath 
to a distance z — —I — —2\ l ~ n . For the constant ion-neutral collision frequency, 
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Figure 4. The electrostatic potential (dashed blue line) and ion fluid flow (solid black 
line) throughout the presheath using the modified mobility limited flow model (left 
axis) . Also shown on the right axis is the estimated ion to electron density ratio Hi / n e 
(dotted red line). 



equation (64) yields e0/T £ 



arccosh(l — z/l) and 



1 - 21/ z) 



(65) 



The profiles for electrostatic potential and ion flow speed through the presheath are 
shown in figure |4| 

Using v = Vi, the stable plasma contribution to the ion-ion collision length is 
plotted in figure [5} We use the criterion that the plasma is collisional if /X l ~ n < 1. 
Figure [5] shows that if one considers only the stable plasma contribution to the ion-ion 
collision length, ions from the entrance to the presheath until about half-way through 
it are collisional. Thus, the ion distribution function in this entrance region should be 
Maxwellian, which agrees with the measurements in Claire et al [33]. However, as the 
sheath is approached, the ion flow speed increases significantly and the ion-ion collision 
length, which is oc v A , becomes much longer than the ion-neutral collision length. This 
suggests that ions in the mid presheath to the sheath edge are nearly collisionless. Thus, 
the stable plasma contribution to ion-ion collisions alone cannot explain the finding by 
Claire et al that the ion distribution function becomes Maxwellian near the sheath edge. 
This suggests that some other mechanism for ion-ion collisions must be present. 

The apparently anomalous ion-ion scattering near the sheath edge in Claire et 
al [33] can be explained by the presence of ion-acoustic instabilities. The instability- 
enhanced collision frequency can be estimated from equation (63) by inputting an ion 
flow speed and ion-to-electron density ratio through the presheath. We use equation (65 ) 
for the ion flow model, which determines k c through equation (A. 10). Estimating the 
density ratio is a bit more difficult because the model in equation (64) takes n.j = n e 
at the sheath edge, but has n, ^ n e into the plasma. Physically, of course, the bulk 
plasma is neutral and neutrality is broken as the sheath is approached. We take for 
an estimate ni/n e w 1/ exp{e[0 — <$>{z = — l))/T e }, which assumes that the ion density 
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Figure 5. Ion- ion collision frequencies normalized to the ion- neutral collision 
frequency through the presheath using the experimental parameters of Claire et al [33] . 
Shown are the stable plasma contribution AJ^ (dash-dotted blue line) , the instability- 
enhanced contribution AJ^* (dashed red line) and the total collision length (solid black 
line). Ions are considered collisional when the total collision length is less than 1 on 
this scale. 



is approximately constant and that the electron density has a Boltzmann drop. The 
resultant density ratio is shown in figure [4] Using these flow and density profiles, the 
instability-enhanced ion-ion collision length is shown in figure [5] 

Figure [5] shows that ion-acoustic instabilities determine the ion-ion collision length 
near the sheath edge, shrinking it by nearly two orders of magnitude. For most of the 
presheath, the ion-acoustic instabilities do not enhance collisions, which are dominated 
by the stable plasma rate. The result is a presheath with three regions: (a) The entrance 
region is collisional because of conventional stable plasma collisions. As the ion flow 
speed increases, these become less significant because oc v A . (b) Ions then become 
collisionless for a region in the mid presheath. (c) Closer to the sheath edge, ion- 
acoustic instabilities enhance ion-ion collisions and ions once again become collisional. 
Thus, the predictions of this model for the plasma parameters of [33] are that ions are 
Maxwellian in both the plasma-facing entrance (a) and sheath-facing exit (c) regions of 
the presheath. In the mid-presheath (b), ions should have a distribution characteristic of 
the collisionless models from section (4.2.1 ). Each of these three regions can be identified 
in the measurements of Claire et al [33]. The predicted ion distribution functions are 
shown schematically in figure [6] A proof that instability-enhanced collisions cause ions 
(and electrons) to evolve to a unique Maxwellian distribution is shown in reference [40J . 
This property holds as long as /u R j "C 1. Here 1j/u R j ~ m e /Mj = 1.3 x 1CT 5 . 

Ion-acoustic instabilities not only enhance ion-ion collisions, but also electron- 
electron collisions. Section 14.2.21 discussed that if electrons are collisionless in the 
presheath, an otherwise Maxwellian distribution function will be truncated beyond a 
certain velocity corresponding to electrons that traverse the sheath and are lost from the 
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Figure 6. Sketch of the predicted ion velocity distribution function at different 
locations in the presheath. This shows (a) a stationary Maxwellian distribution, 
(b) a three-region distribution predicted by collisionless models of the presheath (see 
section 4.2.1), and (c) a flowing Maxwellian distribution. Distribution (a) is expected 



in the plasma-facing region of the presheath where stable plasma ion-ion collisions are 
frequent, (b) is expected in a mid-presheath region where ions are nearly collisionless, 
and (c) is expected near the sheath-presheath boundary where ion-acoustic instabilities 
enhance ion-ion collisions causing them to be collisional again. These correspond to 
the expected distributions in regions (a) , (b) and (c) of figure [5j 



plasma. This truncated Maxwellian has the form of equation (37). However, if A e_e < /, 
electron-electron collisions become frequent enough that they cause this distribution to 
become Maxwellian again. Because momentum is conserved in these electron-electron 
collisions, the resultant Maxwellian will have a slight flow shift towards the sheath. The 
stable plasma contribution to the electron-electron collision length for thermal particles 
is X^ e (v = VTe) = 6.0 m [note that \ l ~ l {v = c s ) = \ e ~ e {v = i>Te)]- This is much longer 
than the presheath length. The stable plasma collision frequency alone leads to the 
prediction that electrons are collisionless in the presheath and should have a truncated 
Maxwellian distribution. 

Claire et al [33] does not present measurements of the electron distribution function 
in their discharge. However, others have measured the electron distribution in similar 
discharges and found it to be Maxwellian, even beyond the energy at which depletion 
should be expected. The first, and most famous, of these measurements was taken by 
Langmuir in 1925 [49J. This apparently anomalous electron-electron scattering near 
the plasma boundaries was later named "Langmuir's paradox" (50]. We recently showed 
that Langmuir's measurement can be understood through the mechanism of instability- 
enhanced collisions from ion-acoustic instabilities [51]. The same physics applies to the 
discharge parameters of Claire et al . Noting that Aj^fv = c s ) = \\^ e {y = t>r e ), the 
relevant ion-ion and electron-electron collision lengths are the same at the sheath edge. 
Thus, ion-acoustic instabilities enhance electron-electron collisions near the presheath 
as well as ion-ion collisions. These are frequent enough to expect a Maxwellian electron 
distribution function close to the sheath [51 J. It is also noteworthy that in the sheath 
itself Ui/n e becomes quite large, which exponentially enhances the instability-enhanced 
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collisions in equation (63). The sheath itself, although spatially narrow, is also collisional 



for ion-ion and electron-electron collisions. 



6. Conclusions 



The central conclusion of this work is that the conventional kinetic Bohm criterion of 
equation Q places undue importance on the low-velocity portion of the ion and electron 
distribution functions. Theoretical literature has generally accepted equation (TP) along 
with the consequence that slow ions dominate the Bohm criterion |22j. This equation 
has been used despite the fact that whenever fi{v z = 0) = 0, or df e /dv z \ Vz=0 = at the 
sheath edge, a divergent integral arrises that renders equation ([!]) unusable. In section |4j 
we considered several common example plasmas in which this occurs. 

We showed in section 2.3 that the primary misstep in derivations of equation ([T| is 



the neglect of a collision operator term of the form J d 3 v C(f s )/v z . Although collisions 
are typically negligible, C(f s )\ Vz=0 ^ when f s (v z = 0) ^ 0, or df s /dv z \ Vz=0 ^ 0. Thus, 
the neglected v~ l moment of the collision operator actually diverges for precisely the 
same distribution functions that cause equation Q to diverge. This term is necessary 
in order to make sense of the divergences that arise in the conventional kinetic Bohm 
criterion. It was also shown that the ion term of equation Q is based on a misapplication 
of integration-by-parts if a general ion distribution function is to be allowed. 

An alternative kinetic Bohm criterion based upon positive-exponent velocity 
moments of the plasma kinetic equation was derived in section |3j The result is 
similar to conventional fluid theories, but where the fluid variables are defined in terms 
of positive-exponent velocity moments of the distribution functions. Developing the 
theory in this manner avoids the possibility of divergent integrals. In contrast to 
the conventional kinetic Bohm criterion, slow ions and electrons have no particular 
significance in this new kinetic formulation of the Bohm criterion. This model was 
compared with equation for several example distribution functions that are common 
in various laboratory plasmas. Equation ([!]) contained divergent integrals for some of 
these examples, while the new theory provided the condition that the ion flow speed be 
supersonic at the sheath edge. 

In section [5j we considered the role of Coulomb collisions in determining the 
ion distribution function in the presheath. Typically ions are assumed to be either 
collisionless or to collide only with neutrals in the presheath because the ion-ion collision 
length in a stable plasma is usually much longer than the presheath length near the 
sheath edge. However, the presheath of plasmas in which Tj/T e <C 1 can be unstable to 
ion-acoustic instabilities. Section [5] shows that these ion-acoustic instabilities shorten 
the ion-ion collision length to such a degree that ions can be considered collisional (to 
ion-ion collisions) near the presheath-sheath boundary and through the sheath. When 
this happens, the ions have a flowing Maxwellian distribution function at the sheath 
edge. A flowing Maxwellian distribution is one example where the conventional kinetic 
Bohm criterion of equation ([I]) diverges. The kinetic Bohm criterion of section [3] simply 



Kinetic Bohm Criterion 



25 



reduces to the fluid-like result of a supersonic ion flow speed for this case. 
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Appendix A. Ion-acoustic instability-enhanced collisions 



Since the ion- acoustic instabilities of equation (62) satisfy ij/ojRj <C 1, the instability- 
enhanced collisional kernel can be approximated by (ID] 



Q 



s—s 
IE 



d 3 k 



kk S[k ■ (v - v')] 5(u>Rj - k • v) exp^t) 



k A ~fj\de(k,u)/du\ 
Corrections to this approximation are of order 0{^f 2 /uj 2 R j) • 



(a.i; 



equation (61), we note that 

de 



m e /Mi <C 1. From 
(A.2) 



Applying the dispersion relation of equation (62), Uj « k-Vj — kc 



n e (l + k 2 \ 2 De ), 



to equation (A.2) gives 
de 2 



du 



n, 



[l + k 2 \ 2 T 



De) 



k 2 C 2 Tli 



A- 4 A 4 
K A De 



(A.3) 



The second delta function in equation (A.I) can be estimated from the more 
elementary form written as a Lorentzian 



<J(wBj-k-v) 



7j 



1 n P 



it (u Rij - k • v) 2 + 7 2 



71 Hi k 2 C 2 



Putting equations (A.3) and (A. 4) into (A.I) yields 

l? s 4 



Q 



s—s 
IE 



2 m. 



k 4 \ 4 



De 



^*-(v-vO] (1 + ( , 2AL)3 . 



,27* 



(A.4) 



(A.5) 



Next, we evaluate 2jt for the convective ion- acoustic waves. Reference [3H] shows 



that the exp(27t) term in equation (A.5) must be calculated in the rest frame of the 



unstable mode. Since the ion-acoustic instability is convective, 



2-ft 



dx! ■ 



v s7 



(A.6) 



'x (k) l v gl 

in which v g = du^/dk is the group velocity, x G (k) is the location in space where 
wavevector k becomes unstable, and the integral dx' is taken along the path of the 



mode. In principle, the spatial integral in equation (A.6) requires integrating the profile 
of 7 and v g , which change through the presheath due to variations in the ion fluid 
speed and the electron density. It also requires knowing the spatial location x G (k) at 
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which each wavevector k becomes excited. In estimating equation (A. 6), we assume that 



changes from spatial variations are weak, and we account for x G (k) by only integrating 
over the unstable k for each spatial location x. Following these approximations we 
obtain 

2Z 7 



2-ft 



(A.7) 



in which Z is a shifted coordinate (with respect to z) that takes as its origin the location 
where the first instability onset occurs. In this case, Z = will be the presheath-plasma 
boundary. The group speed is approximately the phase speed (v g ~ URj/k) for the 
ion-acoustic waves, so 

Z 



2 7 t 



T[m P n; 



kX 



De 



2M t n e X De (l + k 2 \ 2 De ) 3 / 2 ' 



(A.8) 



Returning to evaluating equation (A. 5), we use spherical polar coordinates for k, 
and take the parallel direction along v — v', so that k • (v — v') = k\\(v — v') ~ k\\v. 
Evaluating the fen integral, the k 2 terms become fe 2 = k\. After the azimuthal integral, 
the only nonvanishing components of the tensor Q^ s are the xx and yy components. 
Both components have the same magnitude, which we take as the scalar Ql^ s - We also 
apply the variable substitution k = k^_\r, e . With these, our scalar estimate for the 
collisional kernel becomes 



Q 



Kg; 



IE 



din- 



ts 



exp 



7rm P n,- 



K 



\ De ]/ 2M i n e (l + K^/ 2 



(A.9) 



m s v J Kc (1 + k 2 ^ 3 

in which we have set the lower limit of integration to k c so that only the unstable 
fe are integrated over. The limit k c can be determined from the instability criterion 
Vi — c a f y/l + k 2 > 0, which gives 

V^Pf - 1 > forV^<c s 
, for Vi > c, 



(A.10) 



The k integral in equation (|A.9|) will be approximated as follows. The integrand 

1. Expanding the 



is peaked about the maximum of k 3 /(1 + k 2 ) 3 , which occurs at k 
argument of the exponential about this point yields 



(1 



K 2)3/2 



V2 
4 



V2 



K-l) 



3^2 



(« 



1) 2 + 



(A.11; 



=14 8 32 

Keeping only the lowest order term of this series, we use the approximation k/(1 + 
k, 2 ) 3 ^ 2 w v^2/4 in the exponential. The integrand is then algebraic, and can be evaluated 
analytically 



dn- 



k 3 1 1 + 2k 2 c 

1 + K 2 ) 3 ~ 4(1 + K 2 ) 2 ' 



(A.12) 



With this approximation to the k integral, equation (A.9) becomes 
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